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ABSTRACT 



Context. It is unknown how far dust growth can proceed by coagulation. Obstacles to collisional growth are the fragmentation 
and bouncing barriers. However, in all previous simulations of the dust-size evolution, only the mean collision velocity has been 
0^ ■ considered, neglecting that a small but possibly important fraction of the collisions will occur at both much lower and higher velocities. 

Aims. We study the effect of the probability distribution of impact velocities on the collisional dust growth barriers. 
Methods. We assume a Maxwellian velocity distribution for colliding particles to determine the fraction of sticking, bouncing, and 
fragmentation, and implement this in a dust-size evolution code. We also calculate the probability of growing through the barriers and 
[ the growth timescale in these regimes. 

Results. We find that the collisional growth barriers are not as sharp as previously thought. With the existence of low-velocity 
collisions, a small fraction of the particles manage to grow to masses orders of magnitude above the main population. 
Conclusions. A particle velocity distribution softens the fragmentation barrier and removes the bouncing barrier. It broadens the size 
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distribution in a natural way, allowing the largest particles to become the first seeds that initiate sweep-up growth towards planetesimal 
sizes. 
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C^> 1 . Introduction growth of the larger particle can occur even at very high veloc- 

C*~) ' . ... ,. . , . , ities. In this scenario, the growth initially stalls at the bouncing 

■ Primary accretion is the earliest stage of planet formation, where barrieJ . feut if a ^ number of ^ M j >seed> ides 

• ; tiny, micrometer-sized dust grains in the protoplanetary disk ^ introduced; they rapidly sweep up the smaller particles and 

g . grow into planetesimals of several kilometers in size. In the clas- tQ , ^ Thg h barriers ^ ^ ^ ^ 

^ ■ sic scenario, this happens by the way of incremental growth, in ^ number of , ides and therefore reduce the number of 

^ , which sticking collisions lead to successively arger aggregates destmctive collisions among them . Exactly how the first seeds 

„ . (Weidenschilling 1997; Dullemond & ^Dormrnk 2005). As lab- ^ formed ^ howeyer m ^ ^ 
^ oratory experiments (Blum & Wurm 2008; Guttler et al. 2010) 

and numerical simulations (Zsom et al. 2010; Birnstiel et al. A11 P rior dust coagulation models have until now relied on 
2010) have improved our understanding of the collision process the mean value to describe the velocity at which a collision oc- 
1 and the dust evolution, it has become evident that incremental curs - In realitv > the relative velocities between the particles that 
growth to form planetesimals cannot continue unhindered. arlses because of Browman motion and turbulence does not take 
As particles grow, they decouple more and more from the a sin S le value ' but is better represented by a probability distribu- 
sun-ounding gas, which increases their relative velocities. At tion ' owin g t0 geometrical and stochastic effects. Here, we focus 
the fragmentation barrier, collision energies are high enough to on turbulence since it is the dominating source of relative veloc- 
cause particle destruction, halting the dust growth at centimeter ltv between the small grains below the fragmentation barrier, 
to meter sizes (Braueret al. 2008). Zsom et al. (2010) also in- A general formula for the probability distribution function 
traduced the bouncing barrier, which stops the growth at even (PDF) of particle relative velocities is however unavailable, de- 
smaller sizes. In this case, the collision energies are too low to spite the efforts of many numerical and experimental works, 
cause any particle destruction, but also too high for sticking, with There is tentative evidence that the PDF for particles with large 
growth-neutral bouncing events as the result. Stokes numbers (St ~ 1, those that couple to the driving scales 
Windmark et al. (2012) suggested a sweep-up scenario of the turbulence) is Maxwellian or close to it (Carballido et al. 
where the fragmentation barrier can be circumvented. They 2010; Dittrich et al., in prep.). However, at smaller sizes (where 
found that even though collisions between equal-sized particles particles couple to the Kolmogorov scale) the PDF may be bet- 
generally lead to fragmentation, if the mass-ratio is large enough, ter characterized by wide, exponential tails (Wang et al. 2000; 
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Table 1: Disk model parameters used in the local simulations. 



Parameter Symbol Value Unit 



Distance to star 


r 


l 


AU 


Gas surface density 


E„ 


1700 


g cm -2 


Dust surface density 


Id 


17 


g cm" 2 


Gas temperature 


T 


280 


K 


Sound speed 


Cs 


1.0- 10 5 


cm s -1 


Turbulence parameter 


a 


io- 4 




Solid density of dust grains 




1 


g CITT 3 



Pan & Padoan 2010; Hubbard 2012) Future numerical and ana- 
lytical modelling is desired to refine and interpret these data. In 
this work, as a first step, we assume that turbulent velocities are 
Maxwellian distributed. 

A velocity distribution allows some collisions to result in 
sticking where the average outcome would produce a bounc- 
ing or fragmentation event. This causes the barriers to blur out, 
and might allow for some lucky particles to just by sheer chance 
repeatedly experience only low-velocity collisions and grow to 
larger sizes than the main population. 

In this letter, we show the effect of such a velocity distribu- 
tion in a local dust-size evolution code, not only as a method for 
creating lucky larger particles, but also to see how it affects the 
general dust population. 



With bouncing included (where both particles involved are 
kept unchanged), the sticking efficiency is written as 

f 1 if v < v b 
Ps = \ ifv>v b ' (2) 

where we take a bouncing threshold velocity v b = 5 cm/s, which 
is the upper limit to the bouncing threshold velocity found by 
Weidling et al. (2012). 

In the collision model SF, we only include sticking and frag- 
mentation, and take p s = 1 — pf . In the collision model SBF, we 
account for all three effects, and write the bouncing probability 
as py, = 1 - p % - pf. 

Finally, in the collision model SBF+MT, we include a sim- 
ple prescription for the mass transfer events discussed in detail 
in Windmark et al. (2012). These occur when the particle mass 
ratio is so high that the largest particle can avoid destruction, and 
only the smaller particle is fragmented with a fraction of its mass 
being added to the surface of the larger. To mimic this effect, we 
assume that all fragmenting collisions above a critical mass ratio 
mi/m.2 > 50 result instead in sticking where 10% of the mass of 
the smaller particle is deposited onto the larger one. 

2.2. The velocity distribution 

We assume that all collisions follow a Maxwellian distribution 
characterized by the root-mean-square velocity v rms 



2. Method 

We implement a Maxwellian velocity distribution together with 
three simple collisions models into a local version of the dust- 
size evolution code of Birnstiel et al. (2010). We consider rel- 
ative velocities arising only from Brownian motion and turbu- 
lence, which are the dominant sources for small particles, and 
for the disk properties, we take the minimum mass solar neb- 
ula (Weidenschilling 1977). It is now well-known known that 
the gas and dust mass of the solar nebula might have been much 
higher (Desch 2007), which would increase the coagulation rates 
and the particle size at the fragmentation barrier (Birnstiel et al. 
2011). We refer to Table 1 for all the parameters used in this 
work. 



2.1. Collision models 

In the simple models presented here, the possible outcome of a 
collision is one of the three collision types of sticking, fragmen- 
tation, and bouncing, so that p s + pf + p b = 1, where p is the 
probability of each collision type. In the first model, SF (stick- 
ing + fragmentation), we study the effect on the fragmentation 
barrier, in the second, SBF (sticking + bouncing + fragmenta- 
tion), we study the bouncing barrier, and in the third, SBF+MT 
(sticking + bouncing + fragmentation + mass transfer), we show 
a scenario where growth breakthrough occurs. 

For the fragmentation, we assume that destruction of both 
particles always occurs above a given collision velocity 



Pf 



if v < Vf 

1 if v > Vf 



(1) 



where we take Vf = 100 cm/s as the fragmentation threshold 
velocity as found by Blum & Munch (1993). The fragments are 
put in a size-distribution described by n(m)dm oc mr x ii dm. 



P(.V I Vrms) = \l r— exp -- — 

71 Vfms \ 2v Ls 



(3) 



For the turbulent relative velocity, we use the closed-form ex- 
pressions derived by Ormel & Cuzzi (2007), which for small 
particles (St < 1, here corresponding to m < 10 6 g) can be 
approximately written as v rms = a/9/2 • a St c s , where St is 
the Stokes number and the rest of the parameters are given in 
Table. 1. 

Because of the shape of the velocity distribution, there is al- 
ways a non-zero, if small, chance of sticking during every colli- 
sion event. With the collision prescription discussed in the pre- 
vious section, the total sticking and fragmentation probabilities 
can be written as 



<Ps>(Vrm S )= fV(v)dv 

Jo 

<Pf>(Vrms) = I P(y) dv 
Ji'r 



(4) 
(5) 



In Fig. 1, we plot for the three collision models the integrated 
sticking, bouncing, and fragmentation probabilities as a function 
of particle mass. In the SF model, fragmentation occurs much 
earlier if a velocity distribution is included, but sticking is also a 
possibility at much higher masses. In the SBF model, sticking is 
a possible outcome even orders of magnitude above the bounc- 
ing threshold, but decreases to very low probabilities. At a mass 
of m — 1 g, the sticking probability is 10~ 3 , but fragmentation is 
also rare, so that the relative ratio of the two is of the same order 
of magnitude. This means that growth can still proceed, albeit on 
longer timescales. Finally, in the SBF+MT model, the situation 
is identical to the SBF panel, except that the added effect of mass 
transfer means that collisions above the fragmentation threshold 
can still lead to growth, provided that the mass ratio between the 
particles is large enough. 
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Fig. 1: The probabilities of sticking (green), bouncing (yellow), 
and fragmentation (red) for the three collision models as a func- 
tion of mass. The red-green dashed line in the SBF+MT panel 
represents mass transfer, where fragmentation can turn into mass 
gain for large mass ratios. The dotted lines represent the thresh- 
old masses without a velocity distribution included. 
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Fig. 2: Snapshot of the size-distributions for SF (upper), SBF 
(middle), and SBF+MT (lower) collision models taken after 
t — 5 ■ 10 4 years, with (red) and without (black) a velocity dis- 
tribution. The gray diagonal lines correspond to a total of 1 and 
10 6 particles within a 0.1 AU annulus, respectively. 



3. Results 

In Fig. 2, the dust size distributions of the simulations are given 
at t — 5 ■ 10 4 years for the three collision models, with and with- 
out a velocity distribution. 

3.1. The fragmentation barrier 

In the SF model, dust growth stops in both cases. When there 
is no velocity distribution, this point occurs abruptly, with no 
way of growing larger particles, after they have reached a size 
corresponding to v rms = Vf . 

If a probability distribution is considered, growth is both pos- 
itively and negatively affected. The main peak of the distribu- 
tion shifts to lower sizes, because collisions from the high-speed 
tail of the distribution already start their destructive work before 
the particles reach the barrier (see Fig. 2 inset). However, there 
are also lucky particles that successively experience low-velocity 
collisions, even beyond the nominal barrier. This leads to the tail 
in the size distribution beyond the barrier seen in Fig. 2. 

The probability of a particle reaching a mass m before being 
destroyed can be approximated by assuming that the particles 
with masses around the peak mass, m pea k, dominate the interac- 
tions with the larger particles. A particle therefore has to undergo 
k consecutive sticking collisions in order to grow to a mass of 
m — k ■ m pea k. The survival chance for particles growing from 
"Jpeak to k ■ mp ea k can be written as a product of the sequence 

k 

^survival = f] Ps(m) ■ (6) 
i=l 



In Fig. 3, we plot the cumulative survival probability under dif- 
ferent assumptions. The dashed and solid lines were calculated 
assuming constant sticking probabilities, and the solid line as- 
sumed a velocity/mass-dependent p s . At a mass of m — 5Qm pl .^, 
the relative velocity had increased by a factor of 2, and the stick- 
ing probability had decreased from 0.5 to 0. 1 compared to the 
peak population. 

These numbers compare well to the large-particle tail in 
Fig. 2. The largest particles of masses m - 60m pea k have a den- 
sity decrease of roughly 25 orders of magnitude relative to the 
peak mass, which is roughly the survival probability that we cal- 
culated in the toy-model. Growing to these masses is extremely 
unlikely, but the sheer number of particles ensures that some 
lucky particles make it. 

3.2. The bouncing barrier 

In the SBF model, we note two differences between the two 
cases. One is the discrepancy in the number of small particles. 
Without a velocity distribution, the small particles are depleted 
as there is no mechanism that can create them once the main pop- 
ulation has grown to the bouncing barrier. This depletion is rapid 
at first, but gets less effective when the number of small particles 
drops and the frequency of sticking collisions decreases. If a ve- 
locity distribution is included, collisions in the high-velocity tail 
cause fragmentation and replenish the population of small parti- 
cles. 

The peak of the distribution shifts significantly towards 
larger sizes. The reason for this can be seen in Fig. 1, where 
the low-velocity collisions lead to sticking, but even the highest 
collision velocities are too low to cause any fragmentation. At 
around m = 1 g, the growth halts when the fragmentation prob- 
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Fig. 3: The cumulative survival probability of a particle cross- 
ing through the fragmentation barrier. The dotted and dashed 
lines are calculated by constant p s regardless of particle mass, 
the solid line by a mass-dependent p s from Eq. 4. 



ability increases rapidly, while the sticking collisions are rare. 
This causes a steeper tail of large particles compared to the case 
of the pure fragmentation barrier. 

We find that bouncing collisions never completely halt the 
dust growth, as there will always be a small chance of stick- 
ing. However, the growth timescale may become so large that 
the growth is effectively halted. If we follow an individual dust 
grain, we can write its growth timescale as 

m m 

^growth — j 77" = i \l) 

dm/dt CT Vrms pp p s 

where cr = n(a + a p ) 2 is the collisional cross-section and p p = 
1.4 ■ 1CT 11 g crrT 3 is the midplane mass density of particles that 
it can collide with. The growth will therefore slow down by a 
factor proportional to the decrease in p s relative to unhindered 
coagulation. In the bouncing barrier, this will cause an increase 
in the growth timescale by a factor of 10 3 . Taking p s from Eq. 4 
and the relative velocity prescription of the previous section, we 
find that it takes ~ 10, 000 years for particles to grow to m = 1 
g. If p s were to decrease further, for example owing to a lower 
bouncing-velocity threshold, this timescale would increase cor- 
respondingly. 

3.3. Breaking through the barriers 

In the SBF+MT model, we finally implemented the physics that 
makes it possible for growth also at high velocities. This relies on 
a mass difference between the particles in the disk, but without 
a velocity distribution, such a mass difference never occurs. 

With the velocity distribution included, the bouncing barrier 
can be overcome (see model SBF), and the fragmentation barrier 
is smoothed out (see model SF), which means that a very small 
(10 6 particles in an 0.1 AU annulus) but important fraction of 
particles manage to grow large enough purely by chance. This 
triggers the growth of these few lucky particles by sweeping up 
the smaller grains trapped below the fragmentation barrier. 

Even a single fragmenting collision between two lucky par- 
ticles will create a myriad of fragments that will also be able to 
sweep up the particles trapped below the barrier. This means that 
the rare fragmenting collisions will effectively multiply the num- 
ber of large particles, and with time, even a handful of lucky par- 
ticles can by themselves create a significant population of plan- 
etesimals. 



For sweep-up to occur in the simulations, a very high dynam- 
ical range is required. Though the break-through occurs for such 
a tiny fraction of the population, the sweep-up growth causes a 
rapid increase in both the number and total mass of the larger 
grains. 

4. Discussion and conclusions 

We have found that the collisional growth barriers for dust grains 
are smoothed out and can even be overcome by virtue of a 
probability distribution of relative velocities among dust grains. 
Although improbable, sticky, low-velocity collisions can occur 
at sizes where the mean collisional velocity would lead to only 
bouncing or fragmentation. 

To grow through the fragmentation barrier, a particle needs 
to be lucky and experience low-velocity collisions many times in 
a row, which causes a tail of larger particles to extend from the 
peak of the mass distribution. Assuming a Maxwellian velocity 
distribution, the luckiest particles can grow to around 50 times 
more massive than they would otherwise be. 

The bouncing barrier is even more affected by the existence 
of a velocity distribution, and particles can grow to more than 
three orders of magnitude higher in mass, with the peak being 
shifted by two orders of magnitude. This occurs because low- 
velocity collisions lead to sticking, but even the higher veloci- 
ties are low enough to only cause bouncing. This means that the 
growth can continue unimpededly until the average relative ve- 
locities have increased enough for the fragmentation barrier to 
start to become important. The bouncing barrier is therefore not 
a solid barrier at all, unless the growth timescale becomes too 
long because of the low sticking probability. 

We have found that the low-velocity tail allows some 
lucky particles to grow beyond the bouncing and fragmenta- 
tion barriers, to become the first seeds in the sweep-up sce- 
nario introduced by Windmark et al. (2012). When the effect 
of fragmentation-with-mass-transferis included, these seeds can 
sweep up the smaller particles trapped by the growth barriers, 
and start their growth towards planetesimal sizes. 

Acknowledgements. We would like to thank the referee for helpful com- 
ments. F.W. was funded by the Deutsche Forschungsgemeinschaft within the 
Forschergruppe 759 "The Formation of Planets: The Critical First Growth 
Phase". For C.W.O. support for this work was provided by NASA through 
Hubble Fellowship grant #HST-HF-5 1294.01-A awarded by the Space Telescope 
Science Institute, operated by the Association of Universities for Research in 
Astronomy, Inc., for NASA, under contract NAS 5-26555. 

References 

Birnstiel, T., Dullemond, C. P., & Brauer, F. 2010, A&A, 513, 79 
Birnstiel, T., Ormel, C. W., & Dullemond, C. P. 201 1, A&A, 525, 1 1 
Blum, J. & Munch, M. 1993, Icarus, 106, 151 
Blum, J. & Wurm, G. 2008, ARA&A, 46, 21 
Brauer, F, Dullemond, C. P., & Henning, T. 2008, A&A, 480, 859 
Carballido, A., Cuzzi, J. N., & Hogan, R. C. 2010, MNRAS, 405, 2339 
Desch, S. J. 2007, ApJ, 671, 878 
Dullemond, C. P. & Dominik, C. 2005, A&A, 434, 971 
Guttler, C, Blum, J., Ormel, C. W., & Dullemond, C. P. 2010, A&A, 513, 56 
Hubbard, A. 2012, arXiv:1207.5365 
Ormel, C. W. & Cuzzi, J. N. 2007, A&A, 466, 413 
Pan, L. & Padoan, P. 2010, Journal of Fluid Mechanics, 661, 73 
Wang, L.-R, Wexler, A. S., & Zhou, Y. 2000, Journal of Fluid Mechanics, 415, 
117 

Weidenschilling, S. J. 1977, Astrophys. & Space Sci., 51, 153 
Weidenschilling, S. J. 1997, Icarus, 127, 290 
Weidling, R., Guttler, C, & Blum, J. 2012, Icarus, 218, 688 
Windmark, F, Birnstiel, T., Guttler, C, et al. 2012, A&A, 540, 73 
Zsom, A., Ormel, C. W., Guttler, C, Blum, J., & Dullemond, C. P. 2010, A&A, 
513, 57 



4 



